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I. INTRODUCTION: BASIC FORMULAE 

O ; 

As it is well known, Newtonian gravitational potential, of homogeneous body with constant density p, at point 
, (X,Y,Z) is defined as triple integral over the body's volume: 
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u(X, Y, Z, x, y, z) = [(a; - X) 2 + (y - Y) 2 + (z - Z) 2 ] 1,2 ; (2) 



here G stands for Newtonian constant of gravitation. In spite of almost 400-year-long attempts since Isaac Newton' 
times, the integral (|j) is known in closed form (not in the series!) in quite a few cases 0: 
a) a piece of straight line, b) a sphere, c) an ellipsoid; 
CN| ' note that both cases a) and b) may be considered as particular cases of c). As to serial solution, the integral (1) is 
expressed in terms of the various kinds of series for external points outside the minimal sphere containing the whole 
body (non-necessary homogeneous), as well as for inner points close to the origin of co-ordinates. However in this 
note we do not touch the problem of serial solution and are only interested in exact analytical solution of (|l]). 
' Recently Kondrat'ev and Antonov Q have obtained the rather sofisticate analytical formulas for the gravitational 
potential (and the gravitational energy) of some axial-symmetric figures, namely homogeneous lenses with spherical 
surfaces of different radii. In forthcoming paper |^] we present some new solutions for homogeneous bodies of revolu- 
tion. Here we present the solution of Eq. ([[]) for homogeneous right (=rectangular) parallelepiped, (RP). The great 
simplification available in this case is that boundaries of integration over each of variables x,y,z are fixed and thus 
independent on other two variables B. This allows to get results in the closed form in terms of elementary functions. 
. To this end we widely used Mathematica Q . 

Let RP with center at the origin of Cartesian co-ordinates have lengths of sides 2a, 2b, 2c, so that inside the RP, 
S_h ' —a < x < a, —b < y < b, — c < z < c . Then Eq. (1) may be rewritten in the following symmetrical form: 



r x P cup r z P 

U(X, Y, Z) = Gp I dx dy u(x, y, z)dz, (3) 

J xm J ym J zm 



u{x,y,z) = (x 2 + y 2 + z 2 ) 1/2 ; (4) 

xm = —a — X, xp — a — X, ym = —b — Y,yp = a — Y, zm = — c — Z, zp = c — Z. (5) 

Note that now the boundaries (Q) of integration in Eq. (|3|) depend on X, Y, Z, but this is not so " dangerous" as the 
case of dependence of boundaries on other variables of integration. 

Now we are ready to start calculation of integral (^|) in analytical form. Certainly final result should not depend on 
the order of integration [3] and this may be used to check the calculation. 
First, an integration over z gives: 

uz{X,Y,Z,x,y) = u(x,y,z)dz =[ln(z+l/u)]lZ% l . (6) 

J zm 
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Here u is the same function as in Eq. (4). 
Second, an integration over y gives: 



uyz(X, Y, Z, x) 



yzu 



xarctan h yln(z + 1/u) + zln(y + X/u) 



y=yp 

y=ym 



(7) 



And integration over x gives the final expression for the gravitational potential of homogeneous RP: 



U(X,Y,Z) = Gp 











z—zp - 


y=yp - 


Gp 






v(x, y, z) + v(y, z, x) + v(z, x, y) 


z—zm- 


y=ym- 
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v&, y\ z>) = y'z' W + ^ + V' 2 + Z ' 2 )' 2 ] ~ %r arctan ^ (9) 

x (x + y + z ' ) 

Here co-ordinates X, Y, Z are as given in (5). In Appendix we present the final result (8,9,5) in explicit form. 
These formulas (|[9,5) may look cumbersome and unpractical while in fact they allow to easily (and exactly!) calculate 
gravitational potential of homogeneous RP at arbitrary point inside as well outside the body. The payment for 
this universality does not seem very high. Also we should note that obtaining the result (||,9,5) is non-trivial in 
the sense that direct using of Mathematica's command for integration of (1,2), Integrate [f[x],x], (and especially 
command for definite integral, Integrate[f[x],{x,xm,xp}]), gives enormous complex expressions even leaving some 
integrals uncalculated! 

Note that a physically evident symmetry relative to all three sign changes, X — > —X, Y — » —Y, and Z — > — Z may 

be used to check final formulas. 

Let us consider now the various particular cases. 



II. A PIECE OF STRAIGHT LINE 



This is one of a few known classical results. From Eq.(^), the 3D-potential, at arbitrary point (A, Y, Z), of a piece 
of straight line, PSL, with constant linear density p and length 2c, with center at origin of co-ordinates , PSL laying 
along the Z-axis, is as follows: 



Ui ine (X,Y,Z) = Gpln 



Z+(X 2 + Y 2 + (c-Z) 2 ) 



2^/2 



Z+(X 2 + Y 2 + (c+Z) 2 ) 1/2 



(10) 



The PSL case may be considered as RP with two infinitesimal dimensions dx, dy. There is a scaling invariance: if we 
express all co-ordinates in units of c then we get the universal law coinciding with the potential of PSL with length 
equal to 2. 



III. RECTANGLE 

From Eq.(Q) we may get the 3D-potential, at arbitrary point (A, Y, Z), of a rectangular with constant surface 
density p, with sides of lengths 2c and 2b along the Z and A axes respectively, with center at origin of co-ordinates 
A, Y, Z. Equipotential surfaces are symmetrical relative to the rectangle's plane. The "rectangle" may be considered 
as RP with one infinitesimal dimension dx. 



IV. CUBE 



From all RPs with three dimensions of RP being finite, the case of cube is of the utmost interest. We consider this 
case in detail. First, the gravitational potential of homogeneous cube is symmetric relative to all three co-ordinates 
which may have only even powers. Second, there is a scaling invariance: if we express all co-ordinates in units of cube 
edge half-length, a, then we get the universal law coinciding with the gravitational potential of homogeneous cube 
with edge length equal to 2. 
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A. Some particular points 



Here are values of gravitational potential of homogeneous cube, with density p and with edge length 2a, at some 
particular points (note that here all potentials are given in units of a? G p): 
at the center: 

U(0, 0, 0) = 24 In 1 - 2tt = 9.52017; (11) 



at the vertex: 



at the center of the face: 



at the middle of the edge: 



[/(l,l,l) = 121n^±^-7r = if/(0,0,0); (12) 



157 12 7 2 

U(0, 0, 1) = 4 arctanh — J - - 4 arctan -J-; (13) 



4 1 44 

U(0, 1,1) = 4 In 10 - arctan - - 8 arctan - = 4 In 10 - arctan yyy - n. (14) 

Though the solution we discuss is exact however it is rather complex and is difficult to use, so in the next sections we 
present some serial expansions. 



B. Neighborhood of the center 



For the homogeneous cube with density p and length of edge 2a, the gravitational potential near the center is 
sphcricall-symmetric and quadratic in co-ordinates (note that all co-ordinates are in this paragraph expressed in units 
of a). The deviation from spherical symmetry appears in the terms of fourth and higher orders: 



U(X, Y, Z) = U S = a 2 Gp 



{24 In - 2tt} - {|tt(X 2 + Y 2 + Z 2 )}+ 



{-^(X 4 + y 4 + Z 4 ) + ^(X 2 Y 2 + X 2 Z 2 + Y 2 Z 2 )} + {-^(X 6 + Y e + Z^ + 



3^3 



162^3 v 



(15) 



5 (X 2 Y* + X 2 Z^ + X 4 Y 2 + X 4 Z 2 + Y 2 Z 4 + Y^Z 2 ) - -^X 2 Y 2 Z 2 } 



108 V3 



Here inner braces separate the terms of zeroth, second, fourth and sixth order respectively. Note that no terms with 
odd powers may occur in serial expansion of the potential of cube and more generally the potential of RP. 



C. Potential at main diagonal of cube 

The serial expansion of the gravitational potential at the main diagonals where abs(X) =abs(Y) =abs(Z) = ry/3, 
and r = \J (X 2 + Y 2 + Z 2 ) is distance from the center) of the homogeneous cube up to r 20 is as follows: 



W diag {r) = a 2 Gp Yt'f c 2t I £ J ; c - 24 In - 2tt, c 2 = -2tt, c 4 



3^3' 



_ 136 _ _ 104 _ 54392 . 136 (16) 



° B 27V3 ,CS 567V^' Cl0 2187v^' Cl2 1082565 %/3 ' ° 1A ' I'iDSvT 

_ 842168 _ 8632 

16 r:eom qn^ . > 18 a 70001:0 . Hi ' 20 



55801305 a/3' 4782969 v^P m 1363146165 V3' 

This series represent the actual potential at the main diagonal very accurately; at the "final" point, at the vertex, at 
r = V3, the difference between series and exact solution (4.76016) is only —0.00612. 
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D. Potential at co-ordinate axis 



Taking Y = 0, Z = we get potential at X axis which connects the centers of opposite sides. Even in this case 
exact expression is rather complex, so we present the serial expansion around the center (at point X = 0). Series up 
to X 20 is as follows: 

U{X, 0, 0) = a 2 G p J2Zl° ^X 2 * 5 co = 24 In - 2n,c 2 = - ^ , c 4 = - ^= , 

„ _ 1 „ = 17 13 _ 4307 _ 17 fi7\ 

6 162 V5' 8 1701 71' 10 104976 73' 12 6495390 V5 ' 14 5668704 75' V ' 

_ 174431 _ 229 346867 

16 oc7cifinfi,in , /^T 1 1" w^cfi^mej . /q ' 2U 



2678462640 ^' 7346640384 \/5 ' U 43620677280^' 



E. Approximations by sphere 

There are many ways of comparing the cube, with edge length 2a and density p, with a homogeneous sphere. 

1. First naive approximation is by the "inscribed" sphere of radius R\ = a and with the same density p. This gives 
the next rough estimation for the gravitational energy of cube (here Mi = An/?>pR\ is the mass of sphere): 

„. 3 G M 2 16 2 „ 5 
5 Hi 15 

Note that gravitational potential energy W of any body is of negative sign but in this paper we loosely write all W's 
with positive sign. 

2. Second approximation is by the "inscribed" sphere with the mass equal to cube's mass (density of this sphere will 
be larger than p, density of cube). Then we have the next, less rough, estimation for the gravitational energy of cube 
(here M 2 = 8pa 3 is the mass of sphere): 

w 3GM 2 192_ 2 5 

W 2 - = —Gp 2 a 5 . 

5 a 5 

3. Next approximation is by the sphere whose volume equal to cube's volume. Radius R3 of this "cquivolumc" sphere 
is i?3 = (6/7r) 1 / 3 a. Homogeneous sphere with this radius and with the density p equal to density of cube gives the 
next approximation for the gravitational energy of cube (here M3 = 8pa 3 is the mass of sphere) : 

W 3 = 3 -^l = ^(J)V3 Gp2a 5 = 30.95 G(? a 5. (18) 
R3 5 6 

Gravitational potential at the center of such sphere, 2tt G pR 2 = 27r(6/n) 2 / 3 Gpa 2 — 9.672 Gpa 2 differs from the exact 
value of central potential of cube (see Eq. (pd])) by less than 1.6%. We deduce from this that W3 ([l8]) should be 
rather good approximation to gravitational energy of homogeneous cube. 

4. Radius R4 of sphere with the same density and the same central potential as ones of the cube is [U(0, 0, 0) /2tt] 1 / 2 a, 
with [t/(0,0,0) given by (|TT|) . This gives another good estimation of gravitational energy of cube: 

_ ZGMl = IS (U W ,)\^ = 
5 R 4 15 V 27r J 

To calculate the potential energy of homogeneous cube we used Mathematica to integrate numerically the formula 
(||) over the volume of the cube and obtained the numerical value 

W cube = 30.117 Gp 2 a 5 . (20) 

We note that (|l8|) and ( JXSj ) give very accurate bounds to "exact" value (|2C|). 
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V. POTENTIAL AT THE CENTER OF RP 



We write down the potential in the center of RP: 

U RP (0, 0, 0) = 4 * {abln + 6c In f±f + cd\n f±g- 
a 2 arctan ^% — b 2 arctan r§ — c 2 arctan 

act Oct ca J 

Here (i = (a 2 + fr 2 + c 2 ) 1 / 2 is the main diagonal of RP. Even potential in the center of RP could not be scaled by 
values of a, b, c. Three particular cases are of the larger interest: 
a) cube corresponding to case c = b = a, see (pi]); 

h)long thin stick with square cross-section corresponding to case a » b = c: 

U sUck (0, 0, 0) = 6 2 (-8 \n(b/a) + 12 - 2 tt + 4 In 2); (22) 
c) thin square plate corresponding to the case a << b = c: 

U p i a te{0, 0, 0) = 16 a b hx(V2 + 1) - 2 tt a 2 . (23) 
These formulae will be used further for comparing RP with ellipsoid. 



VI. COMPARISON WITH ELLIPSOID 



It is interesting to compare the homogeneous RP and the homogeneous triaxial ellipsoid with semi-axes A, B, C , 
and with density and central value of gravitational potential as ones of RP. 

The potential in the center of homogeneous triaxial ellipsoid with semi-axes A, B, C and density p is Q : 

f°° ds 

U eU (0,0,0) = 7rpGABC == (24) 

Jo V {A 2 + s) (B 2 + s)(C 2 + s) 

The potential energy of triaxial ellipsoid is || : 

W e u = TTiG M 2 r ds (25) 

10 Jo ^/(A 2 +s)(B 2 + s )(C 2 + S ) 

Here M = 4/3ir p A B C is the ellipsoid's mass. 

From (|4j) and (|2^) , we note the interesting relation between the gravitational potential at the center of homogeneous 
ellipsoid and the gravitational potential energy of the ellipsoid: 

W e u = lu e ii(0,0,0)M eU (26) 
5 

for any semi-axes! 

For the cube we have from (O) and (|20|) the relation 



W cube = 0.39544 U cube (0, 0, 0) M cube , (27) 

which is very close to ( |2^). 

As integral in (|24|) and (25) in general case is expressed only in terms of elliptic integrals, to simplify calculations we 
assume A > B = C and also a > b = c, then we have for the potential energy of ellipsoid of revolution || : 

3 GM 2 , A 

W e ii = f arccosh — . (28) 

5^ 2 ~B 2 ) B K ' 



and for potential at the center: 



(0,0,0) = -^—W eU . (29) 



Here M = 4/3ir pAB 2 is the mass of the ellipsoid of revolution. 

There is no simple relation between [7(0,0,0) and W for RP with arbitrary values of edge lengths and we may use 
relation (26) for rather precise estimation of potential energy of RP with known U(0, 0, 0) from (pij). 
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Potential at vertex of RP 



Potential of RP at vertex is: U (A, B,C) — 1/2 {7(0, 0, 0) that is exactly half of the potential at center of RP! 

VII. POTENTIAL ENERGY OF RP 

To get the potential energy of RP, WRP, we need to calculate triple integral over the volume of body: 

WRP = \p f dX { dY [ U(X, Y, Z, a, b, c)dZ. (30) 



2' 



b 



After some lengthy interactive session with Mathematica, we get the potential energy of homogeneous RP with density 
p and edge lengths 2a, 2b, 2c, which we write down in the following concise form: 

^\f( r . h „\ _1_ flh „ „\ _L fir, „ h\\. fir, h „\ — „-rfi _1_ „.„4 _1_ ~„„3 , ___2. 



WRP = Gp 2 [f(a, b, c) + f(b, c, a) + /(c, a, 6)]; /(a, b, c) = c 5 a 5 + c 4 a 4 + c 3 a J + c 2 a 2 ; 

-c)(d 
a dl 



C5 = 32 . C4 = 32 (d _ dl _ d3) _ l|6 m gd#*±6) _ I|£ l n (^cKd+c) . 



c 3 = - ^ arctan S3 i d = Va 2 + b 2 + c2; df = V« 2 + o 2 ; d3 = Va 2 +c 2 ; (31) 
c 2 = ^(dl - d) + ^(d3 - d) - 16 be 2 ln - 16 6 2 c ln {[=§. 

A. Potential energy of cube 
From (|3l"l), taking c = a, b = a, we get the potential energy of homogeneous cube with edge length 2a: 

WC = 32Gp 2 a 5 { 2 ^~/^~ 1 + J + ln[(V2 - 1)(2 - V3)]} - 30.117Gp 2 a 5 . (32) 

B. Potential energy of thin long stick 

We consider a case a >> b = c that is a case of thin long stick with square cross-section. Leading term in expansion 
of WRP ( |3l| ) gives the potential energy of thin long stick: 

W 8Uck = ^Gp 2 ab A \n^. (33) 
3 o 



(34) 



From this and fl22|) we get for stick: 

Wgtick 



8ao 2 {7 stlcfe (0,0,0) 2' 



C. Potential energy of thin square plate 



Taking one of RP dimension infinitcsimally small, a — > 0, we get, from Eq.(31), the potential energy of thin 
rectangular plate. Note that first non-zero term in expansion is quadratic in a. If we additionally take b — c, then we 
get the potential energy of thin square plate: 

WS = 64Gp 2 b 3 a 2 (ln(%/2 + 1) - — ~ 1 | = 47.5714 Gp 2 b 3 a 2 . (35) 



From ( |23D and (35) we have another limit for relation WUM: 

WS 1 V2-1 



\ab 2 U pl (0,0,0) 2 6 1n(V2 + l) 



.421673. (36) 
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D. Relation between potential energy, gravitational potential and mass of RP 



General picture of relation between potential energy, gravitational potential and mass of RP is shown in the Fig.l. 
We shortly denote it as WU M which means: 

"potential energy ^/(potential at the center U(0) x mass of RP M)". 
For homogeneous ellipsoid WUM = 2/5, see (|26|). 




FIG. 1. WUM for RP with square cross section 2fcx2fc and length 2 a. Abscissas are values of b/a and ordinates are values 
of relation WUM = W/U(0) M that is ratio of gravitational potential energy of RP to product of gravitational potential at 
the center of RP by mass of homogeneous RP. WUM has minimum for cube (b/a = 1), tends to 1/2 at b/a —+ (thin long 
stick) and to | — - ^^^js = .421673 at b/a — > oo (thin square plate). For homogeneous ellipsoid WUM = 2/5. 
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APPENDIX 



Here we present a full explicit expression for gravitational potential of homogeneous rectangular parallelepiped with 
density p and with lengths of edges a, b, c along X, Y and Z axes respectively; X, Y, Z are Cartesian co-ordinates 
of point at which the gravitational potential is calculated, and the origin of co-ordinates coincides with the center of 
parallelepiped. This expression is valid for any point X, Y, Z inside as well as outside the body: 
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U(X,Y,Z,a,b,c)= 

{(-c-Z)-2*ArcTan[((-a-X)*(-b-Y))/(Sqrt[(-a-X)-2+(-b-Y)-2+(-c-Z)-2]*(-c-Z))]- 

(-c-Z)-2*ArcTan[((a-X)*(-b-Y))/(Sqrt[(a-X)-2+(-b-Y)-2+(-c-Z)-2]*(-c-Z))]- 

(-c- Z)-2*ArcTan[((-a-X)*(b-Y))/(Sqrt[(-a-X)-2+(b-Y)-2+(-c-Z)-2]*(-c-Z))]+ 

(-c-Z) ~2*ArcTan [ ( (a-X) * (b-Y) ) / (Sqrt [(a-X) ~2+ (b-Y) ~2+(-c-Z) ~2] * (-c-Z) ) ] + 

(-b-Y) ~2*ArcTan [ ( (-a-X) * (-c-Z) ) / ( (-b-Y) *Sqrt [ (-a-X) "2+ (-b-Y) "2+ (-c-Z) ~2] ) ] + 

(-a-X) ~2*ArcTan [ ( (-b-Y) * (-c-Z) ) / ( (-a-X) *Sqrt [ (-a-X) "2+ (-b-Y) "2+ (-c-Z) ~2] ) ] - 

(-b-Y)-2*ArcTan[((a-X)*(-c-Z))/((-b-Y)*Sqrt[(a-X)-2+(-b-Y)-2+(-c-Z)-2])]- 

(a -X) ~2*ArcTan[( (-b-Y) * (-c-Z) )/( (a-X)*Sqrt [(a -X) ~2+(-b-Y) "2+ (-c-Z) ~2] ] - 

(b-Y)~2*ArcTan[((-a -X) * (-c-Z) ) / ( (b-Y) *Sqrt [(-a-X) "2+ (b-Y) "2+ (-c-Z) "2] )] - 

(-a-X) ~2*ArcTan [ ( (b-Y) * (-c-Z) ) / ( (-a-X) *Sqrt [ (-a-X) ~2+ (b-Y) ~2+ (-c-Z) ~2] ) ] + 

(b-Y) ~2*ArcTan [ ( (a-X) * (-c-Z) ) / ( (b-Y) *Sqrt [ (a-X) "2+ (b-Y) "2+ (-c-Z) ~2] ) ] + 

(a-X) ~2*ArcTan [ ( (b-Y) * (-c-Z) ) / ( (a-X) *Sqrt [ (a-X) "2+ (b-Y) "2+ (-c-Z) ~2] ) ] - 

(c-Z) ~2*ArcTan[( (-a-X)* (-b-Y)) /(Sqrt [(-a-X) ~2+(-b-Y)-2+(c-Z)~2]* (c-Z) )]+ 

(c-Z)-2*ArcTan[((a-X)*(-b-Y))/(Sqrt[(a-X)-2+(-b-Y)-2+(c-Z)-2]*(c-Z))]+ 

(c-Z) ~2*ArcTan [ ( (-a-X) * (b-Y) ) / (Sqrt [(-a-X) ~2+(b-Y) "2+ (c-Z) ~2] * (c-Z) ) ] - 

(c-Z) ~2*ArcTan[( (a-X)* (b-Y)) /(Sqrt [(a-X) ~2+(b-Y) ~2+(c-Z) ~2] *(c-Z) )] - 

(-b-Y) ~2*ArcTan [ ( (-a-X) * (c-Z) ) / ( (-b-Y) *Sqrt [ (-a-X) "2+ (-b-Y) ~2+ (c-Z) ~2] ) ] - 

(-a-X) ~2*ArcTan [ ( (-b-Y) * (c-Z) ) / ( (-a-X) *Sqrt [ (-a-X) "2+ (-b-Y) ~2+ (c-Z) ~2] ) ] + 

(-b-Y) ~2*ArcTan [ ( (a-X) * (c-Z) ) / ( (-b-Y) *Sqrt [ (a-X) "2+ (-b-Y) ~2+(c-Z) ~2] ) ] + 

(a-X) ~2*ArcTan [ ( (-b-Y) * (c-Z) ) / ( (a-X) *Sqrt [ (a-X) "2+ (-b-Y) "2+ (c-Z) ~2] )] + 

(b-Y) ~2*ArcTan [ ( (-a-X) * (c-Z) ) / ( (b-Y) *Sqrt [ (-a-X) ~2+ (b-Y) "2+ (c-Z) ~2] )] + 

(-a-X) ~2*ArcTan [ ( (b-Y) * (c-Z) ) / ( (-a-X) *Sqrt [ (-a-X) ~2+ (b-Y) ~2+(c-Z) ~2] )] - 

(b-Y) ~2*ArcTan [ ( (a-X) * (c-Z) ) / ( (b-Y) *Sqrt [ (a-X) "2+ (b-Y) "2+ (c-Z) ~2] ) ] - 

(a-X) ~2*ArcTan [ ( (b-Y) * (c-Z) ) / ( (a-X) *Sqrt [ (a-X) "2+ (b-Y) "2+ (c-Z) ~2] ) ] }/2- 

(-b-Y) * (-c-Z) *Log [-a-X+Sqrt [ (-a-X) "2+ (-b-Y) "2+ (-c-Z) ~2] ] - 

(-a-X) * (-c-Z) *Log [-b-Y+Sqrt [ (-a-X) "2+ (-b-Y) "2+ (-c-Z) ~2] ] + 

(-b-Y) * (-c-Z) *Log [a-X+Sqrt [ (a-X) "2+ (-b-Y) "2+ (-c-Z) ~2] ] + 

(a-X) * (-c-Z) *Log [-b-Y+Sqrt [ (a-X) "2+ (-b-Y) "2+ (-c-Z) ~2] ] + 

(b-Y) * (-c-Z) *Log [-a-X+Sqrt [ (-a-X) ~2+ (b-Y) "2+ (-c-Z) ~2] ] + 

(-a-X) * (-c-Z) *Log [b-Y+Sqrt [ (-a-X) "2+ (b-Y) "2+ (-c-Z) ~2] ] - 

(b-Y) * (-c-Z) *Log [a-X+Sqrt [ (a-X) ~2+ (b-Y) "2+ (-c-Z) ~2] ] - 

(a-X) * (-c-Z) *Log [b-Y+Sqrt [(a-X) ~2+ (b-Y) "2+ (-c-Z) ~2] ] + 

(-b-Y) * (c-Z) *Log [-a-X+Sqrt [ (-a-X) ~2+(-b-Y) ~2+(c-Z) "2] ] + 

(-a-X) * (c-Z) *Log [-b-Y+Sqrt [ (-a-X) ~2+(-b-Y) ~2+(c-Z) ~2] ] - 

(-b-Y) * (c-Z) *Log [a-X+Sqrt [ (a-X) "2+ (-b-Y) "2+ (c-Z) ~2] ] - 

(a-X) * (c-Z) *Log [-b-Y+Sqrt [(a-X) "2+ (-b-Y) "2+ (c-Z) "2] ] - 

(b-Y) * (c-Z) *Log [-a-X+Sqrt [(-a-X) ~2+ (b-Y) "2+ (c-Z) ~2] ] - 

(-a-X) * (c-Z) *Log [b-Y+Sqrt [(-a-X) ~2+ (b-Y) "2+ (c-Z) "2] ] + 

(b-Y) * (c-Z) *Log [a-X+Sqrt [ (a-X) '2+ (b-Y) '2+ (c-Z) ~2] ] + 

(a-X) * (c-Z) *Log [b-Y+Sqrt [ (a-X) '2+ (b-Y) '2+ (c-Z) ~2] ] - 

(-a-X) * (-b-Y) *Log [-c+Sqrt [ (-a-X) "2+ (-b-Y) "2+ (-c-Z) ~2-Z] + 

(a-X) * (-b-Y) *Log [-c+Sqrt [ (a-X) "2+ (-b-Y) ~2+ (-c-Z) ~2] -Z] + 

(-a-X) * (b-Y) *Log [-c+Sqrt [ (-a-X) "2+ (b-Y) '2+ (-c-Z) "2] -Z] - 

(a-X) * (b-Y) *Log [-c+Sqrt [ (a-X) ~2+ (b-Y) ~2+(-c-Z) ~2] -Z] + 

(-a-X) * (-b-Y) *Log [c+Sqrt [ (-a-X) "2+ (-b-Y) "2+ (c-Z) ~2] -Z] - 

(a-X) * (-b-Y) *Log [c+Sqrt [ (a-X) '2+ (-b-Y) '2+ (c-Z) ~2] -Z] - 

(-a-X) * (b-Y) *Log [c+Sqrt [ (-a-X) '2+ (b-Y) '2+ (c-Z) ~2] -Z] + 

(a-X)*(b-Y)*Log[c+Sqrt[(a-X)-2+(b-Y)-2+(c - Z)~2]-Z]. 

Here the potential is given in units of G p and the Mathematica's language is used except of first " additional" 
lme,U(X, Y, Z, a, b, c) =; Mathematica's designations are corresponding to "usual" mathematical formulas as follows: 
Log[x]= ln(x), Arctan[x]= arctan(x), Sqrt[x]= x 1 / 2 . 
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